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Abstract 

The Beam Energy Scan program has been undertaken at the Relativistic Heavy Ion CoUider 
(RHIC) to search for the QCD critical point. The presence of the critical point is expected to lead 
to non-monotonic behavior of several quantities. Here we report the result of higher moments 
of net-charge distributions for Au+Au collisions at y^Sjviv = 39 GeV as measured by the STAR 
experiment. The STAR results are compared with results from HIJING event generator and 
Hadron Resonance (ERG) Models. 



1. Introduction 

Lattice calculations [1] have predicted the presence of a critical point in the phase 
boundary of hadronic matter with Quark-Gluon Plasma (QGP). The critical point is the 
end point of the first order phase transition. To map the QCD phase diagram and to 
locate the critical point, the Beam Energy Scan program has started in the year 2010 
at the Relativistic Heavy-Ion Collider at Brookhaven National Laboratory. The program 
aims to scan the beam energy from y/sJfN = 5 GeV to 200 GeV corresponding to baryon 
chemical potentials up to as high as 550 MeV and down to 14 MeV[2]. The characteristic 
signature for the CP is the divergence of the susceptibilities of conserved quantities 
like net-charge, net-baryon, net-strangeness [4] and the correlation length(^)[l,3]. These 
quantities are related to the higher moments of the event-by-event distribution of the 
above conserved quantities. The variance, skewness and kurtosis are related to ^^-^ 
and ^^[3] respectively. The ratio of the fourth to the second order susceptibilities are 
related to the product of the kurtosis and variance(iircr^). It is expected that these higher 
moments will be non-monotonic, at the CP, as a function of an experimentally varied 
parameter such as the beam energy [5]. 

We present the preliminary results of higher moments and their product of net-charge 
distribution from Au-|-Au collision at 39 GeV and 200 GeV measured by the STAR 
experiment at the RHIC. 
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2. Analysis Details and Results 



The STAR experiment provides tlie excellent particle identification and large accep- 
tance for event- by- event fluctuation analysis. The Time Projection Chamber(TPC) is the 
main tracking device. It is used to identify the charge particles within full azimuthal and 
±1.8 unit of pseudo-rapidity (77). The event- by- event net-charge is measured at y^J/vw = 
39 GeV for Au -l-Au collisions occurring within 30 cm of the TPC center along the beam 
line. For this analysis, nearly 6 x 10^ minimum bias events are used. The charged particles 
between the range 0.15 < pr < 1.0 GeV/c are measured at \r]\ < 0.5 region. Standard 
STAR track quality cuts are used for this analysis. For the electron and positron iden- 
tification in same pt range, which is a source of the background for the analysis. Time 
Of Flight (TOF) detector is used. For the centrality selection, uncorrected charged parti- 
cles multiplicity measured within \ri\ < 0.5 from the TPC is utilized . To get the average 
number of participant (< Npart >) for each centrality, Glauber model calculation is done. 
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Fig. 1. Net-charge distribution for Au+Au collision at y^sj^ = 39 GcV for 0-5%, 30-40% and 60-70% 
centrality . 

Figure 1 shows the net-charge distribution for centrality 0-5%, 30-40% and 60-70%. 
The net-charge distributions are symmetric around the mean value for different central- 
ities (small skewness values as seen in Fig. 2). The HIJING event generator is used for 
model comparison with data in same 77 and px region. To understand the leptonic con- 
tribution to the higher moments of net-charge distribution, both the model and the data 
based approaches (utilizing the TOF detector) are made. Both these studies show their 
contribution is negligible. 

The mean, standard deviation, skewness and kurtosis of Au+Au collisions at 39 GeV 
are plotted as a function of < Npart > representing different centrality in Fig. 2. These 
moments are fitted with the expected form of the central limit theorem(CLT) to under- 
stand the volume dependence of the higher moments [6]. Both the mean and standard 
deviation increase with < Npart >: where the mean shows linear dependence. The error 
bars include both statistical errors and systematic errors due to the range of multiplicity 
within each centrality bin. Reduction of the centrality bin width effect and the other sys- 
tematics studies are in progress. In more central collisions, the kurtosis shows a deviation 
from the CLT that need more study to understand this effect. 
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To eliminate the volume dependence in moments, the product of the kurtosis and 
variance {Kcr"^) and that of the skewness and standard deviation {Sa) are used. Ka^ 
and Sa are plotted versus < Npart > in Fig. 3 for Au+Au collisions at 39 GeV and 
200 GeV. Here Au+Au 200 GeV results are taken from [7]. The results from the Hadron 
Resonance Gas(HRG) model[8] are consistent with the data within the errors. The results 
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Fig. 2. The mean(top left), standard deviation(top right), skewness(bottom left) and kurtosis (bottom 
right) are plotted versus < Npart >■ The dotted lines are the expected values from the central limit 
theorem. The error bars are statistical and the grey lines represent systematic error due to the centrality 
bin width effect. 
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Fig. 3. The Ka^ (left-panel) and Scr (right panel) of the net-charge distribution are plotted versus 
< Npart > respectively from the data. The filled rectangle represent the HRG model results from [8]. 
The blank rectangles are the HIJING model results for the Au-|-Au at 39 GeV. The red star and the 
blue triangle represent Au+Au 39 and 200 GeV results from the data. The error bars are statistical error 
and black caps represent systematic error due to the centrality bin width effect. 
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from HIJING are also agreement with the data. The Ka^ values are larger than unity 
where the Sa values are close to zero for different centralities. The Ka^ and Sa values 
are similar for Au+Au collisions at ^/snn = 39 and 200 GeV. 

3. Summeiry 

The STAR preliminary results of the higher moments of the net-charge distribution for 
Au+Au collision at ^/snn = 39 and 200 GeV are presented. The results are compared 
with the HRG and HIJING models. Within current experimental uncertainties, the Ka^ 
and Sa values are comparable for collision at 39 and 200 GeV. The results are in agree- 
ment with the HRG model. The results from HIJING model, which does not include CP 
physics, are consistent with those from data. The data analysis and the model studies 
are in progress for the other BES energies. 

References 

[1] R. Gavai and S. Gupta, Phys. Rev. D 78, 114503 (2008). 

[2] J. Clcymans et al., J. Phys. Rev. C 73. 034905 (2006). 

[3] M.A. Stephanov, Phys. Rev. Lett. 102, 032301 (2009). 

[4] M. Cheng et al., phys. Rev. D 79, 074505 (2009). 

[5] M. Stephanov et al., Phys. Rev. D, 60,114028 (1999). 

[6] M.M. Aggarwal et al., PRL 105, 022302 (2010). 

[7] T. K. Nayak for the STAR collaboration, Nucl. Phys. A 830 (2009). 

[8] F. Karsch and K. Redhch, Physics Letters B 695 (2011). 



4 



